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A long time ago, it has been conjectured that a Hamiltonian with a potential of the form x^-\-ivx'', 
V real, has a real spectrum. This conjecture has been generalized to a class of so-called "PT symmetric 
Hamiltonians and some proofs have been given. Here, we show by numerical investigation that the 
divergent perturbation series can be summed efficiently by an order-dependent mapping (ODM) in 
the whole complex plane of the coupling parameter v'^ , and that some information about the location 
of level crossing singularities can be obtained in this way. Furthermore, we discuss to which accuracy 
the strong-coupling limit can be obtained from the initially weak-coupling perturbative expansion, 
by the ODM summation method. The basic idea of the ODM summation method is the notion 
of order-dependent "local" disk of convergence and analytic continuation by an order-dependent 
mapping of the domain of analyticity augmented by the local disk of convergence onto a circle. In 
the limit of vanishing local radius of convergence, which is the limit of high transformation order, 
convergence is demonstrated both by numerical evidence as well as by analytic estimates. 
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I. INTRODUCTION 

The order-dependent mapping (ODM) summation 
method has been introduced in Ref. [1] and initially ap- 
plied to series expansions of an integral virhose perturba- 
tive expansion counts the number of Feynman diagrams 
with four-point vertices, to the quartic anharmonic oscil- 
lator, and to renormalization group (RG) functions in the 
three-dimensional cj)'^ quantum field theory [1, 2]. Later, 
other examples of quantum mechanics and field theory 
type have been studied [3-7]. Some convergence proofs 
have been given in Refs. [8, 9]. However, these exam- 
ples have in common, as a starting point, perturbations 
of Hermitian quantum Hamiltonians. In most of them 
the ODM summation method appears as a higher order 
extension of some variational calculation [1, 10]. There- 
fore, we may wonder whether Hermiticity plays a role 
in explaining some specific convergence properties of the 
method. 

Here, we thus consider a non-Hermitian example, the 
expansion generated by an i addition to the har- 
monic potential. The corresponding Hamiltonian is 
symmetric, that is, symmetric under a simultaneous 
complex conjugation and space reflection. Even though 
this is not quite obvious, such a Hamiltonian has a real 
spectrum, as first conjectured by Bessis and Zinn- Justin 
(1992) and proven in [11, 12]. In this note, wc show 
numerically that the ODM method, suitably adapted to 
the problem, allows a precise determination of the ground 
state energy. The convergence properties also give some 
information about the analyticity as a function of the 
strength of the i perturbation, in particular excluding 
level crossing in some region of the Riemann surface. 



Another potential application of the method is the 
quantum field theory in two and three space dimensions, 
whose renormalization group functions describe the Lee- 
Yang edge singularity of Ising-like statistical models [13, 
14], should longer perturbative series become available. 



II. THE IMAGINARY CUBIC POTENTIAL 

We here consider the spectrum of the simplest example 
of a T^T symmetric Hamiltonian, 



H = - 



1 ^ 

2 d.T2 



1 



(1) 



where 5 is a real positive parameter. Indeed, the Hamil- 
tonian is invariant under simultaneous complex conjuga- 
tion and parity transformation x t-^ —x. In this conven- 
tion, the resonance energies of the cubic potential have 
their branch cut along the negative real g axis, where the 
quantity i ^ becomes real and the particle may tunnel 
through the "cubic wall" of the potential. One then ver- 
ifies that the perturbative expansion of the energy eigen- 
values En for 5 —7^ contains only integer powers of g: 



En{g) = ^En,kg 



k=0 



(2) 



with real coefficients. Moreover, a steepest descent cal- 
culation of the path integral representation of the cor- 
responding quantum partition function [15-17] yields a 
large order behaviour of the form 



En,k 



C {-!)'' k'' A-'' k\ , 



(3) 
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where A — 24/5 is the instanton action. The constant 
C and the half-integer b depend on n. The series is di- 
vergent for all g and when the expansion parameter is 
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not small, a summation of the perturbativc expansion 
is indispensable. Note that the sign oscillation with k 
and some additional considerations already suggest that 
the series is Borel summable (a property proved in [18]) 
and thus that the spectrum, beyond perturbation the- 
ory, is real. This was initially the origin of the conjec- 
ture of Bessis and Zinn- Justin, which was generalized 
to other 'PTsymmetric Hamiltonians [19-21]. More re- 
cently, Pade summability was also rigorously established 
[22]. This proof confirms numerical investigations based 
on the summation of the perturbative expansion by Pade 
approximants [23]. Distributional Borel summability of 
the perturbation series to the complex resonance energies 
for negative coupling g was proved in [24]. 

In this note, we show by a combination of numerical 
and analytic arguments that the ODM method is espe- 
cially well adapted to summing the perturbative series, 
converging for all values of the expansion parameter in 
the complex plane and, in addition, in a region of the sec- 
ond sheet of the Riemann surface. It thus also provides 
additional information about the spectrum for g complex 
on the Riemann surface [25]. 

However, before dealing with the quantum mechanics 
example, as a warm-up exercise, we first consider the 
simple integral 



Let E{g) be an analytic function that has the Taylor 
series expansion 



1 /■"'"°° 

Z{g) = — = / da; exp (-ia;^ - ^^/gx^) 
V 27r J-oo 



(4) 



whose expansion coefficients count the number of Feyn- 

man diagrams that appear in the expansion of the parti- 
tion function tre^^^ corresponding to the Hamiltonian 

(!)• 



FIG. 1. Mapping: example of a function analytic in a cut- 
plane. The cut in the p-plane (l.h.s.) is mapped onto the 
contour V of the A-plane (r.h.s.). 



III. ORDER-DEPENDENT MAPPING 
SUMMATION METHOD 

A. Method 

The ODM summation method [1] is based on some 
a priori knowledge, or educated guess, of the analytic 
properties of the function that is expanded. It applies 
both to convergent and divergent series, although it is 
mainly useful in the latter case [1]. 



E{g) = Y.^,g' 



(5) 



where the equality has to be understood in the sense of a 
(possibly asymptotic) scries expansion. In what follows, 
we consider two functions analytic at least in a sector (as 
in the example of Fig. 1) and only mappings 5 A of 
the form 

^ = P (1 _\)« ' " > 1 ' (6) 

where a has to be chosen in accordance with the analytic 

properties of the function E (a rational number in our 
examples) and p is an adjustable parameter. A general 
discussion of the method can be found in Refs. [1, 2]. 
An important property that singles out such a mapping 
is the following: inverting the mapping, one finds that, 
for g 00, A approaches unity, and has a regular 

expansion in powers of 1 — A (sec the Appendix). 

After the mapping (6), E is given by a Taylor series in 
A of the form 



E[g{\))=Y^P,{p)\'- 



(7) 



where the coefficients Pi{_p) are polynomials of degree I in 
p. Since the result is formally independent of the param- 
eter p, the parameter can be chosen freely. At p fixed, the 
series in A is still divergent, but it has been verified on a 
number of examples (all Borel summable), and proved in 
certain cases [8, 9] that, by adjusting p order by order, 
one can devise a convergent algorithm. The basic idea 
is that p characterizes a "local radius of convergence" of 
the divergent input scries. As the order of the mapping 
is increased, p — > 0. However, the mapping g — g[p, A) 
is constructed so that the circle of convcirgcnce of order- 
dependent radius p is always mapped onto a finite domain 
of the A complex plane. 

At a risk of some oversimplification, we can describe 
the general paradigm of the ODM summation as follows: 
One "pretends" that the divergent input series is in fact a 
convergent series, but with a circle of convergence whose 
radius gradually vanishes as the order of the perturbative 
expansion is increased. Since both p — and also |A| < 
1, the direct expansion in A, with variable p, may then 
lead to a convergent output series even for large g. 

The fc-th approximant E'^'^\g) is constructed in the 
following way: one truncates the expansion (7) at order 
k and chooses p = pk as to cancel the last term, i.e. 
Pk{pk) = 0. Because Pk{p) has k roots (real or complex), 
in general one chooses for p the largest possible root (in 
modulus) pk for which Pj-ip) is small (but see the more 
detailed discussion of Sec. IIIB). 

This leads to a sequence of approximants 
k 

E^''Hg) = J2PMX\9,Pk) with Pk{pk)=0, (8) 
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where A = \{g,p) is obtained by inverting the rela- 
tion (6). In the case of convergent series, it is expected 
that pk has a non- vanishing Hmit for k ^ oo. By con- 
trast, for divergent series it is expected that pk goes to 
zero for large k as 



-i/k 



(9) 



The intuitive idea here is that pk is proportional to a 
"local" radius of convergence. 

The following remarks are in order, (i) Alternatively, 
one can choose the largest roots pk of the polynomials 
P'kip) for which Pk is small. Other mixed criteria involv- 
ing a combination of Pk and Pk can also be used. Indeed, 
the approximant is not very sensitive to the precise value 
of Pk, within errors. Finally, Pk+i{pk)^^^^ , as well as 
the values of Pk{p)X'^ in the neighbourhood of pk, give 
an order of magnitude of the error, (ii) In the ODM 
method, the main task is to determine the sequence of 
Pk- Indeed, once the pk are known, for each value of g 
the calculation reduces to inverting the mapping (6) and 
simply summing the Taylor series in A to the relevant 
order. 



B. Optimal convergence analysis 

We here give a heuristic analysis [1] of the convergence 
of the ODM method that shows how the convergence can 
be optimized. This will justify the choice of the class of 
zeros of the polynomials Pk or and provide a quanti- 
tative analysis of the corresponding convergence. 

In this note, we consider only real functions analytic 
in a cut plane with a cut along the real negative axis (see 
Fig. 1) and a Cauchy representation of the form 



9' -9 



zO) 



(10) 



up to possible subtractions to ensure the convergence of 
the integral for g — oo. For later purpose, we introduce 
the notation 



A(5) = luiE{g + i{)). 



(11) 



Moreover, we assume that the function A asymptotically 
fulfils 

A(<?) (X (-5)-^-V/^ A>Q. (12) 

For the two examples discussed below, this exponential 
asymptotic form can be derived by steepest descent cal- 
culations [15-17]. The function E{g) can be expanded in 
powers of g: 

E{g) =Y,Eig' with E, = - h{g) . 

(13) 



The assumption (12) then leads to a large order be- 
haviour of the form 

Ek oc {-Ayrik + b+l) ^ {-A)-''k''k\ (14) 

k—>oo 

as assumed in Sec. 11. We now introduce the mapping 
(6), 



A 



9 = P 



(1-A)" 



(15) 



The cut on the real negative axis of the gi-plane is then 
mapped onto the contour f (Fig. 1). 

The Cauchy representation then can be written as 

'^iam _ 1 I^y^i9m 



2in Jf A' - A TT Jr 



A- A' 



(16) 

where F is the image of the upper-part of the cut on the 

real negative axis under the mapping A = X{g) at fixed 
p, with a segment of the real A negative axis, 1/(1 — a) < 
A' < 0, and the complex contour Im(f ) > ending at 
A = 1. The contours can be deformed if the function 
E{g) has analyticity properties beyond the first Riemann 
sheet. 

We expand, omitting the dependence on p in A, 



E{g{X))^J2PMH9)]' 



(17) 



where the approximants are obtained by the truncated 
expansion 

k 

E{g{X))^J2PMH9)Y, (18) 
and pk is chosen so that Pk{pk) ~ 0. Moreover, 

Pe{p) = l£^A{giX)). (19) 

For — )• oo, the factor A~*^~^ favours small values of |A|, 
but for too small values of A, the exponential decay of 
A(g(A)) takes over. Thus, the remainder value of the 
polynomial Pk{pk) can be evaluated by the steepest de- 
scent method. The ansatz 



Pk^R/k, R>0, 



(20) 



implies that the saddle-point values of A are independent 
of k and g{X) — t- 0. Thus, A{g) can be replaced by its 
asymptotic form (12) for (/ — > 0_, except for g of order 1 
and thus A close to 1. In what follows, we set 



R 



(21) 



since this is the only parameter (and it is independent of 
the normalization of <?). The behaviour of Pk{pk), with 
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Pk being given by Eq. (20), is then given by the leading 
saddle point contributions: 

Pk{pk)=0{e'"') with a = ^(l-Ar-lnA. (22) 

At a zero of Pk{pk), several leading contributions can- 
cel but each contribution is expected to give an order of 
magnitude of the error. 

The saddle point equation is 

""'^^^ " " + Am] = . 

(23) 

For p large, the equation has a unique solution (for /i — >■ 
0, A — l//i), which is real negative and also exists for 
all values of with 1/(1 — a) < A < 0. Then we note 
that at this saddle point, as a function of /U, 



do; 
d/i 



(1 - A)° 
Am^ 



>o. 



(24) 



since the saddle point value of A is negative. This sug- 
gests decreasing ji as much as possible to improve the con- 
vergence, and thus, taking the zero pk with the smallest 
modulus. The exponential rate a corresponding to the 
saddle point vanishes for 



— (l-A)«-ln(-A)=0, 



(25) 



and this defines a special value jic of the parameter /z. 
Some numerical values, obtained by combining Eqs. (23) 
and (25) are displayed in Table I. For /i > /ic, this con- 
tribution to Pk{pk) grows exponentially while for p < Pc 
it decreases exponentially. 

The further analysis then somewhat depends on 
whether a is smaller or larger than 2. If a > 2, there 
exists a region on F where a is positive and convergence 
is only possible if the initial integration contours F or f 
can be deformed in what corresponds to the second sheet 
of the fimction E{g). More generally, the contribution of 
the contour from the neighbourhood of A = 1 plays an 
important role except, again, if it is possible to deform 
F in the first expression (16) in such a way that A = 1 is 
no longer on the contour but inside it. 

Two cases need to be distinguished: 

(i) One can decrease p until other, complex, saddle 
points or maxima of the modulus on the contour give 
contributions with the same modulus. Relevant zeros 
of the polynomials Pk and P^. then correspond to the 
cancellation between the different saddle point contribu- 
tions. If then p > pc a is positive and the new approx- 
imants eventually also diverge and just yield a better 
behaved asymptotic expansion, li p < pc, cr is negative 
and the contributions to Pk{pk) decrease exponentially 
with k (this implies the possibility of contour deforma- 
tion), and, as we show later, the method converges for 
all values of the parameter g. This also implies that the 
mapping (6) removes all singularities of E{g), obviously 



an exceptional situation but illustrated by the example 
of Sec. IV. 

(ii) Generically, the mapping (6) does not cancel all 
singularities of the function E{g) and if it is possible to 
decrease p up to pc, a vanishes and then other contri- 
butions of order at most exp (Cfci-i/") (in particular 
coming from the integration near A = 1) appear and the 
relevant zeros of Pk and correspond to the cancella- 
tion of these contributions. One can further decrease p 
but the convergence properties are unlikely to improve 
since contributions of order exp (C survive and 

can no longer be reduced. 

Returning to the expansion (17), at g fixed, from the 
behaviour of pk one infers 



1- A 



and so 



A^ 



exp 



1/a 



R 



(26) 



(27) 



If the saddle point with A < is a leading saddle point 
and p 7^ pc [case (i)], convergence is dominated by the 
behaviour of e'^'^ and this irrespective of the value of g. 

If the integral is dominated by contributions of order 1 
in the sense of case (ii) , p must be chosen in the range p < 
Pc to avoid exponential divergence, but then to decrease 
the factor A*^ one should choose R, and thus p, as large 
as possible. This implies choosing p = Pc- 

In a generic situation, we then expect the behaviour of 
the contributions to Pk{pk) to be dominated by a factor 
of the form 



exp(cA;^-i/") 



(28) 



and the domain of convergence depends on the sign of 
the constant C. For C > 0, the domain of convergence is 



|5| <i?C7-«[cos(arg(5)/a)]«. 



(29) 



For a > 2, this domain extcinds beyond the first Riemann 
sheet and requires analyticity of the function E{g) in the 
corresponding domain. 

For C < 0, the domain of convergence is the union of 
the sector |arg(5)| < 'Ka/2 and the domain 



|5|>i?|C|-"[-cos(arg(5)/a)]". 



(30) 



Again for a > 2, this domain extends beyond the first 
Riemann sheet. 



IV. EXAMPLE OF THE ia;=' INTEGRAL 

We now study, as an example, the integral 



1 

V27r J-oo 



(31) 
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TABLE I. Values of /ic and Ac- 



3/2 
4.031233504 
0.2429640300 



2 

4.466846120 
0.2136524524 



5/2 
4.895690188 
0.1896450439 



3 

5.3168634291 
0.1699396648 



4 

6.1359656420 
0.14003129119 



The function has a divergent series expansion in powers 
of g, 

Z{9) = 1 - ^ff + + 0{g') , (32) 

and the coefficient A in (3) is given by the non-trivial 
saddle point 



a; + 5 ^/g =0 => x = 



2i 

71' 



and so the saddle-point value of the integrand is 
e^/'? with A = I . 



(33) 



(34) 



Moreover, the function Z{g) has a convergent large g 
expansion of the form 



Z{g)=g-'/'J2'^9- 



e/3 



where the evaluation of the integral yields 



(35) 



zo = 6^/^ = 1.1212331717419689582 . . . (36) 

3 ^(3) 

This suggests the mapping (see the Appendix) 

9 = (37) 
and the introduction of the function 

^{X,p) = {l-X)-y'Z{g{X)). (38) 

Setting in the integral x = s^l — A, one obtains 

1 / \ 

<f){X,p) = ^= / ds exp (-is^ + iAs^ - i -/Xps^l 
V27r J -00 ^ ' 



(39) 



where 



1 1 /■ dA 



^ 2m Jr A*^+i 

X / ds exp(-is2+ iAs2- l^s^) . (40) 

J —00 



The first terms are 

</.(A,p) = l+Q-^p) A 



3 35 



385 



p^]X^ + OiX") . (41) 



8 48^"^ 1152' 



The asymptotic behaviour of Z for g ^ 00, 

^(5) -^0 5"'/' (42) 



is obtained for A ^ 1, 

^(5)~5-'V/V(l,p)- 

It follows 

cl){l,p) = zop-^'^. 



(43) 
(44) 



We use the following values of pk in order to calculate the 
first ODM approximants (j)'^''^ [1, p) to </>(!, p) according 
to Eq. (8), but with the condition Pl^{pk) = (zeros of 
the derivative) except for fc = 1 where we choose the 
condition Pi(pi) = 0. These read (for fc = 1, 2, 3, 4) 

2.4, 1.0909, 0.7058 ± 0.1866 i 0.5894 ± 0.2633 i . 

(45) 

The approximants to zq defined in Eq. (36) yield (we give 
only the real parts when pk is complex) 



p^'^(f)^^\l,p) = 1.15709373, 



(l,p) = 1.26825944, 



ol/6^(2) 



p^l^(j)^^\l,p) = 1.18358984, 
pi/^^^Hl-P) = 1.08048484. 



(46a) 
(46b) 
(46c) 
(46d) 



Surprisingly, using only a minimum number of terms 
from the input series, a rather good approximation to 
the strong-coupling limit is obtained. 

Figures 2, 3, 4 and 5 show that for Pk{p), in a rescaled 
variable r = kp/A, there is an essentially fc-independent 
range where Pk (p) and (p) are small. This is also the 
region in the neighbourhood of which the zeros of Pk 
and are located. Quite generally, we choose for pk the 
complex zeros of P^(p) with the largest real part in this 
range because they are easy to systematically identify. 

The expected theoretical value, obtained by expressing 
that the real negative and the two complex conjugate 
solutions of the saddle point equation (23), 



A= -0.1897168... , 

A = 0.8448584 . . . ± i 1.386261 . 



(47a) 
(47b) 
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TABLE II. Values of and errors in the ODM approximants 
in the limit g — >■ oo. For the estimated convergence rate of the 
output series obtained by the ODM, {SkY^'^ should approach 
a constant as — >■ cxa. 



k 


5 


10 


15 


20 


25 


30 


k Re(oL ] 
A 


3.3612 


4.3300 


4.3450 


4.4210 


4.4954 4.5365 


fclm(pfcj 
A 


0.7187 








0.4544 


0.2984 


0.2260 




0.4910 


0.6149 


0.6480 


0.6526 


0.4972 


0.6819 


k 


35 


40 


45 


50 


55 


60 


kRe(pk) 
A 


4.54721 


4.5509 


4.5580 


4.5698 


4.5840 


4.5884 


fclm(pfc) 
A 


0.1800 


0.4522 


0.3706 


0.5774 





0.4390 




0.7008 


0.7064 


0.7004 


0.7019 


0.7232 


0.7323 





FIG. 4. The plot shows Pk{p) as a function of r = kp/A for 
k = 60, for the model problem (31), demonstrating that a 
range where Pk{p) is small persists in higher transformations 
orders. 



0,001fr 

0,0002- 




FIG. 5. The plot shows APk{p)/k as a function of r = kp/A 
for k — 60, again for the model problem (31). We demonstrate 
that the favourable region where both Pk{p) and Pk{p) are 
small persists in higher transformations orders. 



FIG. 2. The ODM polynomial Pk{p) as a function of the 
scaled variable r = kp/A for k = 30. The model problem (31) 
is being studied. 




FIG. 3. The plot shows APl.(p)/k as a function of r = kp/A 
for k = 30, for the model problem (31). 



yield contributions that have the same modulus (which 
allows for compensation) is = R/A k, 4.62987613 . . . 
and then the rate of convergence is predicted to be of the 
form 0.775*^ from the corresponding estimate (22). 

The numerical investigation for g — >■ oo extrapolated 



to fc = oo leads to (see Table II and Fig. 6) 
lim A:Re(pfc)/yl = 4.65±0.03, 



(48) 



yielding a value of /i significantly smaller than the value 
5.31 ... of Table I in section III B and consistent with the 
theoretical expectation. 

To calculate the ODM approximants, we fit pk com- 
bining the theoretical asymptotic form and the calcu- 
lated values of pk in the range k < 60. This is achieved 
by adding a slightly ad hoc small correction term to the 
asymptotic formula. 

For 5 — >■ oo, we define the following quantity which 
estimates how well the asymptotic strong-coupling limit 
is approximated by the ODM approximants of order fc, 



ODM,k 



(49) 



Some data are then displayed in Table II and Fig. 7. 
From a fit to the numerical data in Fig. 7, One infers a 
geometric convergence of the form (0.78 ± 0.02)^^ quite 
consistent with the theoretical prediction. When g is 
finite and for fc — > 00, the additional factor coming from 
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A*' has the form (see Sec. Ill B) 



exp - 



52/5 ) 



(50) 



This factor can never cancel the geometric convergence 
factor and, thus, the ODM approximants converge on the 
whole Riemann surface. 



0,030 0,035 0,040 0,045 



FIG. 6. kKc(pk)/A, odd-cvcn averaged, as a function of 1/k, 
for the model problem (31). The plot illustrates how the 
asymptotic limit given in Eq. (48) is approached for k ^ oo. 



0,020 0,025 0,030 0,035 0,040 0,045 



FIG. 7. The plot shows (4/519)^^^''"^''' as a function of 1/k. 
We illustrate the asymptotic rate of convergence of the ODM 
approximants to the strong-coupling limit for the model prob- 
lem (31), in higher transformation orders beyond k = 19. 



V. QUANTUM MECHANICS 

A. Strong coupling limit of the cubic potential 

Wc now consider the spectrum of the Hamiltonian (1), 
which can be obtained by solving the time-independent 
Schrodinger equation 

-^r{x) + (ix2 + 1 7^ x^) V(x) = E ^Pix) (51) 

with appropriate boundary conditions. The path inte- 
gral representation of the corresponding partition func- 
tion has a structure that bears some analogy with the 
integral (31) discussed in Sec. IV. 



The spectrum has an expansion in integer powers of g 
with real coefficients. For example, for the ground state 



384^ ' 



(52) 



Instanton calculus, based on a steepest descent evalu- 
ation of the corresponding path integral, allows us to 
calculate the classical action relevant for the large order 
behaviour (3). One finds 



5 



(53) 



For g large, from the scaling properties of Eq. (51), one 
infers that the ground state energy Eq = Eo{g), as a 
function of g, for large g, has an expansion of the form 



,1/5 



E 



U9 



-21/5 



(54) 



This suggests the mapping (see the Appendix) 
A 



(1 - A)5/2 



i?o = (l-A)-i/20(A,p). (55) 



In what follows, we again first concentrate on the be- 
haviour of i?o(ff) for 5 — >■ oo, which is related to 4>{l,p). 
More precisely, taking the strong coupling limit A ^ 1, 
one finds 



-^(I.P) 



(56) 



We now proceed as for the model problem discussed in 

Sec. IV. Figures 8 and 9 show in a rescaled variable, that 
there is a range, which is /c-independent up to corrections 
of order A;~^/^, where Pk{p) and its derivative are small. 
This is also the region in the neighbourhood of which the 
zeros of Pk and are located. We have first determined 
the complex zeros pk of (p) with the largest real part in 
this range. Again, these zeros have been chosen because 
it is easy to identify them systematically. 

A fit of the numerical data for orders A; < 55 and an 
extrapolation to fc — > oo then yields 



fcRc(Pfc) 
A 



and thus 



5.5 ±0.2 



R 



6.0 ±0.5 



- = 5.5 ±0.2. 
A 



(57) 



(58) 



to be compared with the expected maximum value fj, = 
4.895 ... of Table I. 

The numerical relative rate of convergence is compat- 
ible with an exponential of the form 0.05 x 0.94*^ and 
yields the estimate 



0.3727 ±0.0003. 



(59) 
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FIG. 8. We show Pk{p) as a function of t — kp/A for the 
ODM approximant of order fc = 55 to the ground state energy 
of the cubic Hamiltonian (51). 




FIG. 9. The plot displays API. (p) /fc as a function of the scaled 
variable r = kp/A for the ODM approximant of transforma- 
tion order fc = 55 to the ground state energy of the cubic 
Hamiltonian (51). 



0.22 0,2.1 0,24 Q,25 0.26 0.27 0.2S 



FIG. 10. kRe{pk)/A as a function of 1/k^^^. The range for 
the abscissa is the interval (0.20,0.29), and an extrapolation 
of the data to fc — >■ oo is compatible with Eq. (57). 



However, the value of fi would lead eventually to a di- 
vergence of order 1.2*^. This indicates that the form (57) 
cannot be used asymptotically. If we take into account 
the expected value /i = 4.895 . . . and again fit Re(pfc), we 
find 



^= 4.895690188-^ 



(60) 



FIG. 11. i5o,fc, even-odd averaged, as a function of 1/fc^/^. 
As fc — > oo, the energy converges — under oscillations — to the 
point indicated as 1/k^^^ — >■ 0. The plot illustrates the find- 
ings of Eqs. (61) and (62). 



form, we indeed find a smoother convergence with 



,1/5 



^ 0.37255 ±0.00004, 



(61) 



at order 55 to be compared with the value obtained by 
solving the Schrodinger equation Eq — 0.372545791 [27]. 
The errors can be fitted by the form 

■^{Eo.k + Eo^k-i) — -Eoas. (62) 
- 0.00092 exp[-0.292 fc^/^] cos(4.02 {k^/^ - 1.20)), 

where the odd-even averaging improves the fit of the er- 
ror (and slightly the convergence) by cancelling some sub- 
leading oscillatory correction. The form of the oscillatory 
factor in (62) is just a simple guess (see Fig. 11). 



B. Finite value of the coupling parameter 

For g finite, the additional factor (27) governing the 
convergence of the ODM approximants in transformation 
order k is 



with 



^ ~exp( -57^ 



ss 3.535. 



(63) 



(64) 



As a check, we have summed the series for g = 5. The 
result is 



Eo{5) = 0.60168393320519(6), 



(65) 



and the convergence rate of the order of e"^'^'^ '^^^ . Com- 
paring with the convergence for 5 — ^ 00, one finds a vari- 
ation of the coefficient of k^^^ of 1.81 to be compared 
with 3.5/^2/5 ^ 1.86. 

For g = 1, the result obtained from a numerical solu- 
tion of the Schrodinger equation yields 



a form we use in the remaining calculations. With this 



£;n(l) = 0.53078175930417667113556, 



(66) 
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with all digits being significant. The ODM at order 55 
yields (we give its explicit numerical value for illustra- 
tion) 

Eo{l) = 0.53078175930417667113565 . (67) 

The analysis of the convergence is consistent with the 
theoretical form with a coefficient of k^^^ equal to 3.85 
compared to the theoretical expectation of 3.82. 

Finally, to compare with the Fade summation of 
Ref. [23], we give results for g = 288/49, which corre- 
sponds to = 1/7 in Ref. [23], 

A; = 54 , Eoi^) = 0.6127381063888 , 

A-Bo(^) = 5.524167213054 , (68) 

A; = 55, £;o(lf) = 0.6127381063891 

A£;o(^) = 5.524167213068 . (69) 

where AE^ ^ A9{Eo - At order 54, in Ref. [23] 
one finds AEq 5.5241683 and, at order 192, AEq = 
5.52416721306031. The theoretical formula for the ODM 
methods predicts an error for AE of the order of a few 
times 10~^^ for k = 192, confirming the improved conver- 
gence achievable by use of the ODM. For completeness, 
we have verified that AEq = 5.5241672130602221133 . . . 
by numerical solution of the Schrodinger equation. 

C. The imagiriciry part on the real negative axis 

The ODM method converges on the real negative axis 
and thus the Stieltjes property [22], that is, the positivity 
of the discontinuity, can be verified. For g ^ 0- with 
Im{g) > 0, the imaginary part is known analytically from 
semi-classical considerations (instantons) [27]: 

fi p24/5g / icq ^ \ 

ImEoig) = -^'-^ (l + ^g + Oig'')) , (70) 

where we recall that we are using the cubic Hamiltonian 
here in the convention (1). For [^j — )• oo and lm{g) = 0+, 
one finds 

£;o(.g ^ -oo + iO) - (0.301396 + 10.218977) l.gji/^ (71) 

Above the negative real axis, the antiresonance energy 
has a positive imaginary part. For g = —21.6, at order 55 
we find Eq = 0.5540544 + 0.3513991 to be compared with 
the value Eq = 0.554053519 0.3514017781 given in [27]. 
Similarly, for g = -5 + iO, Eq = 0.4338906 + 0.18385821, 
for g = -1 + iO, Eo = 0.4425200442 0.0155179271. In 
fact, the imaginary part is a simple positive decreasing 
function providing a smooth interpolation between the 
two asymptotic forms (see Fig. 12). 

D. Convergence and singularities 

Combining Eqs. (62) and (63), one finds that the con- 
vergence domain is given by 

Re(£?-2/5) > -0.0826 . (72) 



8- 
7- 

o 

6- 

<> 

5- 

4- * 
3- o 

o 

o 



i-, 1 . 1 . 1 . 1 . 1 . 1 . 1 . r — ' — i^-' — f 

2 4 6 8 10 12 14 16 18 20 

FIG. 12. The function — lri(lmiJo(</)) as a function of —g. 

This contains the sector |argg'| < 57r/4, that is, the whole 

Riemann sheet and a sector of the second sheet. More- 
over, in the neighbourhood of 5 = 00, one finds the ad- 
ditional domain 

\g\ > 510.0 [- cos(2arg5/5)]5/'. (73) 

In the union of these two domains, the function Eo{g) 
is analytic and thus free of level-crossing singularities. 
In particular, the expansion (54) has about 0.0826 as a 
radius of convergence. 

VI. CONCLUSIONS 

In this article, we have explored the properties of 
the order-dependent mapping (ODM) for the summation 
of the divergent series originating from non-Hermitian 
Hamiltonians. In these cases, the resulting resonance en- 
ergies may be complex, and their behaviour may be stud- 
ied as a function of the complex coupling parameter g. In 
Sec. II, we discuss the motivation for the current study 
and the two physical problems treated here, which arc 
the cubic anharmonic oscillator in quantum mechanics, 
and a model problem that counts the number of Feynman 
diagrams that occur in the perturbative expansion of the 
partition function of the cubic anharmonic oscillator. 

In Sec. Ill, we recall some important properties of the 
ODM summation method, and in particular, we describe 
how the saddle-point approximation helps in determin- 
ing asymptotic estimates for the parameter pk that are 
needed in order to construct the ODM in higher orders, 
for the convergence rate of the resulting approximants in 
the strong coupling limit, and for the analyticity domains 
of the approximants as the complex phase of the coupling 
parameter is varied. 

In Sec. IV, we discuss the model problem defined in 
Eq. (IV) which is a simple integral that counts Feynman 
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diagrams in a cubic theory and gives rise to a divergent 
series. We demonstrate that good approximations to the 
strong couphng limit can be obtained on the basis of 
few perturbative terms [see Eq. (46)]. In Sec. V, the 
analysis is first carried over to the symmetric cubic 
potential, with similarly encouraging numerical results, 
before demonstrating the applicability of the method to 
manifestly complex resonance energies (see Sec. VC). 

The basic idea of the ODM is that one pretends that 
the divergent input series, which is to be summed, is 
in fact a convergent series, but with a circle of conver- 
gence whose radius gradually vanishes as the order of 
the perturbative expansion is increased. This algorithm, 
together with the representation of the coupling param- 
eter g given in Eq. (6), leads to a very rapid convergence 
rate of the ODM transforms, in part because the dou- 
ble expansion in p and A implied by Eq. (6) represents 
g as a function of two variables p and A whose modu- 
lus is bounded by unity for approximants constructed in 
the strong coupling limit. If one additionally makes a 
suitable change of variable (and function), in order to 
incorporate the information about the strong coupling 
asymptotic expansion (see the Appendix), then one may 
approximate the strong coupling behaviour using only 
very few perturbative input data. 
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Appendix A: Pcirameters in the ODM 

Let the strong- coupling asymptotic expansion for the 
quantity E{g) be known, 

i^(5) = /E^"5-"/", 5^oo, (Al) 

rt=0 



a property shared by the two examples we have discussed 
here (Eqs. (35), (54)) but also by the quartic anharmonic 
oscillator [1] and all perturbations to the quantum 
harmonic oscillator. We consider the conformal mapping 
g = p\{l — A)~". This transformation maps the real 
positive fif-axis onto the finite A interval [0, 1]. 

For 5 — > 00, A ^ 1 and g has an expansion at A = 1 of 
the form 



5-^/" = ^A„(l-A)"+i (A2) 

n=0 



with Ao = p-i/". The function 



</.(A) = (1 - \r^E{g{\)) (A3) 
has then a Taylor series expansion at A = 0, 



0(A) = ^(/.„A", (A4) 

n=0 



as well as at A = 1, 



<^(A) = ^^„(1-A)" (A5) 

n=0 



with i^f) = eoP , where eo is the coefficient defined in 
(Al). This last property explains, to a large extent, the 
good convergence of the method even for 5 — )• oo. 
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